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This work is an ab initio study of the 2p 3 4 S^ 2 , and 2 Dy 2 s , 2 states of C~ and 2p 2 3 Po,i,2, 
1 D2, and 1 So states of neutral carbon. We use the multi-configuration Hartree-Fock approach, 
focusing on the accuracy of the wave function itself. We obtain all C~ detachment thresholds, 
including correlation effects to about 0.5%. Isotope shifts and hyperfine structures are calculated. 
The achieved accuracy of the latter is of the order of 0.1 MHz. Intra-configuration transition 
probabilities are also estimated. 

PACS numbers: 32.10.Hq, 31.15.aj, 31.15.ac, 32.10Fn, 31.15.ag 



I. INTRODUCTION 



ls 2 2s 2 2p 2 'S 



J=2 



Negative ions have always attracted broad attention 
from the scientific community [U [2] • They challenge both 
the experimentalist and theoreticians. The first because 
they are weakly bound, and therefore fragile, and because 
they do not possess a lot of features allowing measure- 
ments. The latter because the binding of an extra elec- 
tron is granted only by a arrangement of the electrons in 
a highly correlated system |3|. Moreover, the fact that 
the electrons in negative ions are bound by a short range 
potential confer them unique properties. 

C~ is the lightest negative ion to have two bound 
terms: its ground state 4 S° and the 2 D° excited state 
which both arise from the 2p 3 configuration. The level 
diagram of the states studied in this work is given in 
Figure [T] 

Carbon is among the most abundant components in 
the universe and a key element in life chemistry. The 
carbon negative ion is attracting in astrophysics and 
atmosphere physics since nitrogen-like 2p 3 A S° — 2 D° 
forbidden lines are recognized as useful transitions for 
abundances determination 4-6 . It has also recently 
been suggested by Le Padellec et al. [7] that nega- 
tive ion could intervene in astrophysical reactions. The 
photodetachment cross-sections of the C~ have been re- 
peatedly studied, both theoretically [5] and experimen- 
tally [S] , for photon energies addressing valence electrons 
and core electrons jTol [TT] . Recently, an isotope separa- 
tion method was tested by Andersson et al. [12] . based 
on the isotopic dependence of the Doppler shift of the 
C~ detachment thresholds in an accelerator. 

A binding energy of 1.262 119(20) eV for the C~( 4 S°) 
has been measured by Scheer et al. [13] who could not 
improve the old value of 33(1) meV for the C~( 2 D°) de- 
tachment threshold, measured by Feldmann [14] , The 
fine-structure of the 2 D° state is not known. On the the- 
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FIG. 1. Level diagram of the C and C~. The fine structure 
of the 3 P state of the neutral carbon and the 2 D° state of C~ 
are magnified ( x 100) . 



oretical side, very accurate carbon electron affinities were 
obtained with coupled-cluster-based methods [HHIZ]- 

The structure of the C~ has not been studied thor- 
oughly and especially little is known about the 2 D° mul- 
tiplet. In laboratory plasmas, lifetimes of the order of the 
ms were measured for the C~( 2 D°), the electron detach- 
ment being principally caused by the black-body radia- 
tion and, to a lesser extent, to collisions [18]. Significantly 
longer lifetimes could be reached in the cold and diluted 
interstellar media where molecular anions have already 
been detected [15] , However, the C~( 2 D°) is, for various 
reasons, very difficult to study experimentally. In this 
context, a firm theoretical knowledge of this system is 
particularly precious. 

Elements from boron to fluorine are the next targets af- 
ter beryllium in the working line of "exact" calculations. 
High accuracy can be achieved for systems with up to 
four electrons using wave functions expanded in explicitly 
correlated gaussian or in Hylleraas coordinates 
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Although the precision that can be achieved for atoms 
with more electrons is limited by the complexity of the 
electron correlation mathematical treatment, the ground 
states of the second period p— block atoms from B to F 
are satisfactorily described by a non-relativistic approach 
on top of which relativistic corrections arc added. 

A critical benchmark quantity for highly correlated 
models is the isotope shift (IS) on the electron affin- 
ity (54) that is doubly sensitive to correlation effects: 
through the negative ion structure and through the 
specific mass shift parameter. The multi-configuration 
Hartrcc-Fock method has been successfully used for cal- 
culating the IS on the 54 of O [23], S P3 and CI [25]. 

The purpose of the present work is double. Our princi- 
pal objective is to obtain the crucial informations about 
the C - electronic structure for stimulating experimen- 
tal research on the 2 D° state. Therefore, we focus on 
quantities that are especially difficult for experimental- 
ists to measure: isotope shifts, hyperfine structures and 
absolute transition probabilities. As for the energy sep- 
arations themselves, we do not try to compete nor with 
the observation, nor with the previous coupled-cluster 
calculations. We instead use these reliable reference data 
for assessing the quality of our computational procedure. 

Our second objective is to obtain non-relativistic (NR) 
wave functions as accurate as possible using the standard 
tools of the ATSP2K package [55] • For getting the best 
estimation of the accuracy, we choose to use the same sys- 
tematical construction of our C and C~ models, avoiding 
any arbitrary compensation of the "additional" electron 
correlation of the negative ion compared to the neutral 
atom. 

In Section [il] we present large scale numerical multi- 
configuration Hartree-Fock (MCHF) calculations (Sec- 
relativistic calculations estimated using the 
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Breit-Pauli approach (BPCI, Section II F) [37] and rel 
ativistic configuration interaction based on the Dirac 
equation (RCI, Section |IIG[) 28 . In Section [TTT 



we 



present accurate results for hyperfine structures (Sec- 



tion |III A ), total energies including the fine structure 
(Sec tion [n iB]), mass polarization shift parameters (Sec- 



tion 



IIICy~oTC 2p 2 3 P, 1 D, 1 S and C" 2p 3 4 5°, 2 D°. 
In Section |IIID| we present the Ml and E2 transition 
probabilities within the 2p 2 and 2p 3 configurations of C 
and C~. 



II. COMPUTATIONAL METHOD 

A. The MCHF expansion 

The multiconfiguration Hartree-Fock (MCHF) varia- 
tional approach consists in optimizing the one-electron 
functions spanning a configuration space and the mixing 
coefficients of the interacting configuration state func- 



tions (CSF) [27] for describing a given term 

y(~/LSM L M s ) = Ci$(rfiLSM L M s ). (1) 

i 

B. Hyperfine interaction 

The level hyperfine structure is caused by the interac- 
tion of the angular momentum of the electrons (J) and of 
the nucleus (I) , forming the total atomic angular momen- 
tum F = I + J. The theory underlying the computation 
of hyperfine structure using MCHF wave functions can be 
found in references [2"5H3"T] . It is possible to express the 
non relativistic hyperfine interaction in terms of the J- 
independent orbital (a;), spin-dipole (a s d), contact (a c ) 
and electric quadrupole (b q ) electronic hyperfine param- 
eters defined as [5S] 
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N 

b q = (TLSM L M S \ 2C£ ] (i)r~ 3 \TLSM L Ms) , (5) 

and calculated for the magnetic component Ml = L and 
Ms = S [32] • The diagonal hyperfine interaction energy 
correction is usually expressed in terms of the hyperfine 
magnetic dipole (Aj) and electric quadrupole (Bj) con- 
stants as follows 

W(J,J) = Aj£ 

3C(C + 1) -4J(J +!)■/( J + 1) 
+Bj 8/(21 -1)J(2J-1) ' (6) 

The first three parameters ([2|, ([3]), and Q contribute 
to the magnetic dipole hyperfine interaction constant 
through 



Aj = A 1 j + A s / + A c j, 



with |33j 



l J- G »T ai LJ(J + 1) 
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(8) 
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SL{2L- 1)J(J+ 1) 



A c , 



A*/ 
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while the last one (b q ) constitutes the electronic contri- 
bution to the electric quadrupole hyperfine interaction 



Bj= —G q Q b q 

6(L-J) 2 - 3(L ■ J) - 2L(L+1)J(J 
L(2L - 1)(J + l)(2J + 3) 



1) 



(.11) 



Expressing the electronic parameters ai, a sc i and a c in 
atomic units (% 3 ) and hi in nuclear magnetons (^jv)j 
the magnetic dipole hyperfine structure constants Aj 
are calculated in units of frequency (MHz) by using 



G L 



95.41067. Similarly, the electric quadrupole hy- 



perfine structure constants Bj are expressed in MHz 
when adopting atomic units (oq 3 ) for b q: barns for Q and 
G q = 234.96475. The expectation values of the angular 
momenta scalar products are given by 

(L ■ J) = [J(J + 1) + L(L + 1) - 5(5 + l)]/2, (12) 
(S ■ J) = [J(J + 1) - L(L + 1) + 5(5 + l)]/2, (13) 
(S ■ L) = [J( J + 1) - L(L + 1) - 5(5 + l)]/2. (14) 

when calculated with non-relativistic LSJ wave func- 
tions. The expression for the off-diagonal hyperfine in- 
teraction, depending on the hyperfine constants 
Bj j_i and Bj j_ 2 , are developed in reference [3Tj . Hib- 
bert [32] gives the expressions of Aj t j-\ in terms of the 
hyperfine parameters (pn5|)- 



C. The isotope shift 

The first order isotope shift on an energy level is de- 
composed in a field shift (or volume shift) and a mass 
shift [21]. The first is proportional to the change in nu- 
cleus rms radius and change of the modified electron den- 
sity at the origin. It is negligible in our context. 

The energy corrected for the first order mass shift, on 
the other hand, can be estimated using [35] 
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where to is the electron mass, M is the bare nucleus mass, 
Z?oo the infinite mass nucleus and 
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The first term contains the normal mass shift 
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and the second one is the specific mass shift (SMS). The 
mass polarization parameter, S sms , has the dimension of 
an inverse square length. 



D. Transition probabilities 

The Einstein A^f coefficient of spontaneous emission 
is defined as the total probability per unit of time for an 
atom in a given energy level i to make a transition to any 
state of the energy level / [35] • 

A transition between levels of same parity is forbid- 
den in the electric dipole approximation, being in general 
many orders of magnitude lower than an allowed transi- 
tion. Two interactions of the same order of magnitude 
can contribute to the appearance of such transitions: the 
dipole magnetic and the electric quadrupole radiation- 
matter interactions. At the non-relativistic level, a dipole 
magnetic transition (Ml) is governed by the electronic 
magnetic dipole operator that is 
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In the mono-configuration approximation, the above ma- 
trix element is non-zero only between states of the same 
configuration and LS. This selection rule is relaxed by 
configuration and LS mixings, the remaining constraints 
being that Jf = J, , Jj ± 1, and that and "J/ have the 
same parity. For its part, an electric quadrupole (E2) 
transition rate is proportional to the electric quadrupole 
moment matrix element 
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(19) 

the sum running on all spatial electron coordinates k. 
Neglecting the LS term mixing, a necessary condition 
for A~E2 to be non-zero is that 5/ — Si =0, \Lf — Li\ < 2, 
\Lf + Li\ > 2 and that the atomic parity is conserved. 



E. Non-relativistic calculations 

We first select a zero-order set of CSFs, the multi- 
reference (MR) . For all studied states it is the set of sin- 
gle and double excitations of the main configuration to 
the n — 2, 3 shells. All the CSFs interacting to first order 
with the MR are selected and we choose the reverse order 
for the subshell coupling [24]. The orbital active set is 
defined as the set of all orbitals characterized by quan- 
tum numbers n < n max and I < l m ax, and is denoted 
\n max l max ~\. We first performed calculations defined in 
the spaces [4/] to |~12/c], denoted MR-I[rt ma:I .Z ma;I ,] . 

For each active space [lOfc] , [life] and [12fc] , we order 
the configurations according to their weight [37 . We 
then construct several new MRs following this hierarchy, 
independently for each state and active set, by selecting 
the minimum group of configurations that add up to a 
certain percentage p of the total wave function. Those 
multi-references are denoted MR p for each given MR- 
I wave function. Unsurprisingly, the MR p sets are not 
sensitive to the used active set. 
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Configuration-interaction calculations are performed 
on each MR p -I CSF sets, p being limited to 99.8% for 
(T( 4 S ) and to 99.3% for CT ( 2 L>°). An example of the 
convergence of the calculations with the number of corre- 
lation layers (n — 2) and p is given in Figure[2j It presents 
results on the 2p 2 3 P state of neutral carbon. The black 
squares show ^SSsms with 



SS S 



versus 



|S sms (MR-I|>iZl) - S sms (MR-I[12fe])|,(20) 



SE = \E(MR-l\nl]) - E(MR-I\l2k])\ (21) 

for n = 6 — 11. Similarly, the white squares com- 
pare the E and S sms convergences of the MR p -I|~12fc] , 
p = 99.0-99.9, results toward the MR 99 .95-I|~12fc] model. 
The energy always decreases along a sequence of increas- 
ingly large calculations, according to the variational prin- 
ciple, while the S sms of MR p -I \nl] calculations decreases 
with n and increases with p. SS sms and SE show a close 
to linear correlation, i.e. the angular coefficients in the 
log-log figure is ~ 1. The slope of this correlation is 
slightly smaller than one for the convergence in n and is 
about 10 — 20 for the convergence in p, as can be seen 
from the offsets in the log-log figure. Similar behaviors 
where found in the open-core CI calculations of S/S~ [24] 
and Cl/Cl - [25]. In general, we can make the following 
observations [35]: the CSFs that are important for the 
energy are accordingly important for the S sms , and the 
Ssms value is more sensitive to the choice of the CSF 
space than to the orbital basis set. 
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FIG. 2. Bi-logarithmic plot of the convergence of the C( 3 P) 
mass polarization expectation value versus the corresponding 
energy, in Hartrees (Eh). The coordinates of the black squares 
are the differences (in absolute value) between the results of 
the MR-I\nl] and MR-I[12fc] calculations, n = 6 - 11 from 
left to right. Similarly, the coordinates of the white squares 
show the convergence of MR P -I[12fc] toward MRg9.g5-I[12fc] , 
p = (99 - 99.9)% from left to right. 

The energy, S sms and hyperfine parameters calculated 
in this work are presented in Table [T] for neutral carbon 



and in Table 151 for C~. It should be stressed that the dif- 
ferences between the results obtained with p = 99.9 and 
p = 99.95 are close to the expected numerical accuracy. 
We estimate that the dominant error in the neutral car- 
bon calculations is due to the nZ-truncation of the active 
set. This is not the case in the C~ calculations for which 
the p-truncation is the most limiting. 

From Tables [I] and [n| we observe that the differences 
between the MR p -I|~llfc] and MR p -Ifl2fe] results do not 
depend strongly on p. In fact, the error made by re- 
porting the impact of the 12 th shell on the calculation 
with p = 99, on the results of the largest MR P -I[11A;] is 
smaller than 4 10~ 7 Eft on the energy, than 2 10~ 6 a^ 2 on 
S sms and than 9 10~ 5 a^ 3 on the hyperfine parameters. 

Using this observation, we add a correction for the 13 th 
correlation layer and for the I = 8 orbitals. First MCHF 
calculations are performed on the MR-I[13fc] and MR- 
I|"12f| CSF spaces, fixing all one-electron radial functions 
at the MR-I[12fc] level and varying only the new orbitals. 
We use the so optimized orbitals in MRgg-I|~13Z] CI cal- 
culations, omitting the 13/ subshell in the active set and 
using the multi-reference obtained with the [~12fc] active 
set. The two contributions (higher n and I) are of same 
order of magnitude, as far as the energy is concerned. 
Still, the additional correlation layer tends to dominate in 
C~ while the additional angular flexibility has the largest 
impact in the neutral carbon calculations. 

As we already mentioned, the convergences in either n 
or p, of a state energy and S sm s are monotone and cor- 
related (see e.g. Figure [2]). This fact could help strongly 
for extrapolating the energy and Ssms value. However, 
even for a two-electron system, the precise behavior of the 
energy convergence with the principal quantum number 
in the high n-limit is unknown [351 ED] ■ Froese-Fischer 
used the following extrapolation function for studying 
four electrons systems [4Tj 



o 4 /(n - Sn) 4 + a 5 /(n - Sn) 5 + a 6 /(n - Snf (22) 



the 04, 0,5 and ae parameters, and Sn being chosen such 
that a.4 < and \ai\ ~ \a 5 \ ~ |a 6 1 . Fitting then = 10 — 12 
results of Table [I] to equation ( [22] ) for extrapolating to 
n — » oo, we obtain —37.84465 E/, for the energy of the 
C( 3 P) state. This procedure does not extrapolate to I — > 
oo. The error on the extrapolation is of the order of 
10~ 5 Eh and the truncation in I of about 10~ 4 Eft. We 
are in fair agreement with the non-relativistic energy of 
-37.8450 Eft estimated by Chakravorty et al. 02]. To our 
knowledge, the values of Table [I] are the best ab initio 
estimated energies, even without the n = 13 and I — 
8 corrections. Finally, let us mention that Sarsa et al. 
[4"3"] calculated the S sm s expectation value for the carbon 
3 P state using the Monte Carlo (MC) approach with an 
explicitly correlated wave function and obtained S sm s — 
-0.38(2) Og 2 . Our final estimated value (-0.40314 a^ 2 ) 
falls a bit outside the statistical MC error bars. 
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TABLE II. Results of the MCHF and CI calculations performed for the C~ 2p 3 4 S° and 2 D° terms. The energies E are in E h , 
the Ssms in a ( y 2 and the hyperfine parameters in a ( y 3 . The final values are the results of the larger \12k] calculations on which 
the impact of the 13 th shell and I = 8 orbitals have been additively transferred. 



model 


ls 2 2s 2 2p 3 


4go 




Is 2 


2s 2 2p 3 2 D° 












nl p 


E 


S S ms 




E 


S sms 


ai 






b q 






MCHF 










MCHF 








HF 

4 

5 

6 

7 
8k 
9k 
10k 
Ilk 
12k 


-37.708844 

— 37.862042 
-37.876688 
-37.884109 
-37.887227 
-37.888691 
-37.889449 

— 37.889853 
-37.890085 

— 37.890213 


-1.60530 

— 0.56521 
-0.56168 
-0.54166 
-0.54675 
-0.54785 
-0.54912 

— 0.54975 
-0.55017 
—0.55042 


0.0 

0.33257 
0.18050 
0.55357 
0.43564 
0.41389 
0.46136 
0.45976 
0.45517 
0.45714 


-37.642589 

— 37.810185 
-37.827492 
-37.836040 
-37.839993 
-37.841966 
-37.842933 

— 37.843464 
-37.843751 

— 37.843927 


-1.54597 

— 0.51517 
-0.51410 
-0.49077 
-0.49006 
-0.49306 
-0.49444 

— 0.49534 
-0.49583 
—0.49614 


2. 
2. 
2. 
2. 
2. 
2. 
2 
2. 
2 
2. 


.35963 
.27217 
28209 
.25909 
.23484 
22493 
.21796 
.21563 
.21476 
.21462 


0.47193 
0.50749 
0.54068 
0.52516 
0.52138 
0.52273 
0.52136 
0.51911 
0.51910 
0.51922 


0.0 

0.23109 
0.23971 
0.35738 
0.31616 
0.30328 
0.32378 
0.32252 
0.32312 
0.32144 


0.0 

0.10808 
0.12912 
0.15402 
0.18659 
0.19303 
0.19962 
0.20107 
0.20291 
0.20238 






CI 










CI 








Ilk 99.0 
99.3 
99.5 
99.7 
99.8 


-37.890143 
-37.890301 
-37.890474 
-37.890640 
-37.890692 


-0.54957 
-0.54774 
-0.54414 
-0.54204 
-0.54132 


0.46140 
0.45872 
0.44761 
0.43942 
0.43823 


-37.844634 
-37.844971 


-0.48786 
-0.48375 


2. 
2. 


.20369 
.20081 


0.51222 
0.51100 


0.30714 
0.30452 


0.21266 
0.21527 


12k 99.0 
99.3 
99.5 
99.7 
99.8 


-37.890271 
-37.890440 
-37.890603 
-37.890769 
-37.890822 


-0.54982 
-0.54797 
-0.54438 
-0.54227 
-0.54151 


0.46325 
0.46098 
0.44919 
0.44080 
0.43945 


-37.844826 
-37.845165 


-0.48798 
-0.48386 


2. 
2. 


.20269 
.19983 


0.51165 
0.51040 


0.30450 
0.30176 


0.21262 
0.21519 


13; 99.0 


-37.890429 


-0.55003 


0.46636 


-37.845003 


-0.48828 


2. 


.20286 


0.51161 


0.30547 


0.21300 


Final 


-37.890980 


-0.54172 


0.44257 


-37.845343 


-0.48415 


2 


.20000 


0.51036 


0.30273 


0.21557 



TABLE III. Relativistic corrections (fiEh) to the total energies evaluated by BPCI calculations (see text). 

Model C ls 2 2s 2 2p 2 CT ls 2 2s 2 2p 3 

nl p A P„ g A P 2 L D2 go «Sg /2 z Dl/2 z D°,/2 

Main spectroscopic terms only (see text) 

lOfc 99 -14 437.25 -14 362.42 -14 239.95 -14 288.85 -14 256.62 -14 200.89 -14 192.92 

life 99 -14 409.01 -14 334.21 -14 211.76 -14 330.16 -14 296.56 
life 99.5 -14 455.75 -14 380.91 -14 258.35 -14 331.96 -14 299.28 

With additional spectroscopic terms (see text) 



lOfc 
life 



99 
99 



-14 437.25 -14 362.63 -14 240.02 -14 288.99 -14 256.62 -14 129.53 -14 201.14 -14 193.15 
-14 409.01 -14 334.42 -14 211.82 -14 330.30 -14 296.56 
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TABLE IV. Relativistic corrections to the S S ms specific mass shift parameter of each considered state evaluated by MCHF-RCI 
MR-SD calculations. Results are presented in [la^ 2 . 



model 






C 2p 2 








C" 2p 3 






3 Pi 




L D 2 




l no 
^3/2 


^3/2 


-^5/2 










mono-configurational 








RCI 


-47 


-47 


-47 


-49 


-53 


-58 


-61 


-61 


DF 


-876 


-560 


65 


-243 


-233 


-91 


-57 


-40 










multi-configurational, 


mono-reference 








3 


62 


190 


375 


316 


343 


449 


439 


442 


4 


4555 


4769 


5143 


4933 


4859 


5080 


4993 


5011 


5 


3996 


4228 


4636 


4388 


4413 


4595 


4678 


4700 


6 


3837 


4067 


4473 


4259 


4376 


4230 


4131 


4154 


7 


3918 


4157 


4579 


4338 


4324 


4305 


4180 


4210 


8 


3939 


4175 


4593 


4337 


4338 


4329 


4272 


4304 










multi-configurational, multi-reference 








3 


98 


225 


412 


351 


390 


441 


442 


449 


4 


4750 


4967 


5349 


5128 


5134 


5201 


5106 


5128 


5 


4167 


4402 


4817 


4559 


4633 


4713 


4794 


4820 


6 


3999 


4233 


4647 


4426 


4570 


4347 


4241 


4268 


7 


4081 


4323 


4752 


4506 


4509 








final 


4102 


4341 


4766 


4506 


4523 


4446 


4382 


4418 



F. Breit-Pauli calculations 

A first way to include relativistic effects is to use the 
Breit-Pauli Hamiltonian that includes the 1/c 2 relativis- 
tic correction operators to the non-relativistic atomic 
Hamiltonian [27] . 

Since the radiative transitions we consider are essen- 
tially authorized by L and S mixing, we need to have 
a good description of the term mixing. On the other 
hand, it is the calculation of the scalar relativistic effects 
that is needed for estimating the relativistic effects on the 
electron affinity since the fine structures of the involved 
species are usually known experimentally. We therefore 
choose two distinct Breit-Pauli models. The first BPCI 
CSF lists are used for the term separation and detach- 
ment thresholds corrections while the second approach is 
used for the transition probabilities calculations. 

Focusing on the correlation, we merge the MRgg-IflOfc] 
lists of the studied terms for both C~ and C. Then we 
extend this model by adding the CSFs interacting to first 
order with the CSFs 2p3p LS, LS = 3 D, 3 S, 1 P for C 
and the CSFs 2p 2 3p LS, LS = 2 P°, 2 F°, 4 D°, 4 P° 
for C~. For the neutral carbon we test the impact of 
additional correlation on the relativistic corrections by 
using the MR g9 -I|ll&f| and MRgg.s-Ifll/c] spaces. We 
finally diagonalize the Breit-Pauli Hamiltonian in those 
CSF spaces using the corresponding active sets optimized 
in the non-relativistic calculations. The relativistic cor- 
rections to the energy are summarized in Table III We 



see that the effect of the additional LS mixing on the 
energy levels is so small that only the corrections on the 
fine structures are meaningful. 



For the reasons expressed above, we also perform BPCI 
calculations that focus on term mixing. For each LSJ 
and active set \nV\ (nl = 4/ — 12fc for C and 4/ — 8k 
for C - ), the MR 98 -I list is merged with the MR-I set 
obtained using the reference containing all allowed LS 



couplings of the 2s — > 3d, 2p — > 3p and 2s 2 
tations from the main configuration. 



2p 2 exci- 



G. Relativistic configuration interaction 
calculations 

We use essentially the same method as in [33]. First, 
we perform reference MCHF calculations with all single 
and double configuration excitations (SD) of the ground 
state in active sets ranging from [3d] to |~8fc] . The result- 
ing non-relativistic radial orbitals P n i(r) are then con- 
verted to Dirac spinors using the Pauli approximation 

P nK (r) =P n l(r) 

^ , s a f d 
Q nK (r) = - ^ + 

where a is the fine structure constant and n is defined 



Pni{r) 



(23) 
(24) 



-I - 1 
I 



when 
when 



3=1 + 1/2 
j = l- 1/2 



(25) 



We finally perform the corresponding RCI calculations 
using the set of SD excitations of the main configuration. 
Larger configuration sets are explored by means of non- 
relativistic CI and RCI calculations using the references 

MR(C~) = {ls 2 2s 2 2p 3 , ls 2 2s 1 2p 3 3d 1 } , (26) 
MR(C) = {ls 2 2s 2 2p 2 , ls 2 2s 1 2p 2 3d 1 , ls 2 2p 4 }. (27) 



TABLE V. Relativistic corrections on A-^ [kHz//ijv] and B/Q [kHz/barn] of each considered state evaluated by MCHF-RCI 
MR-SD calculations. 



in r\(\ pi 

(n \ 

\i ('max ) 




2o 2 


3p 




2p 2 l D 




2p 3 


Ago 




2v 3 
&p 


2 D° 






Ri 


An 


R. 
-D2 


A 2 


B 2 


A3/2 


B3/2 


An In 


B3/2 


Ar In 


-°5/2 












mono-configurational 










RCI 


93 


2 


118 


90 


-141 


-82 


2 


0.1 


69 


855 


-24 





DF 


-158 


-10 


164 


-177 


-26 


150 


-122 


-0.2 


-464 


550 


371 















multi-configurational, mono- 


reference 










3 


-221 


-51 


-52 


38 


-166 


-130 


-404 


-0.4 


-62 


1317 


-113 


26 


4 


-226 


-31 


157 


-111 


-23 


61 


-301 


-0.1 


-18 


1276 


50 


24 


5 


-256 


-34 


152 


-121 


-13 


70 


-322 


0.0 


-16 


1280 


68 


23 


6 


-254 


-35 


163 


-123 


-3 


64 


-319 


0.2 


-19 


1356 


73 


27 


7 


-279 


-39 


158 


-135 


6 


71 


-357 


-0.1 


-19 


1393 


77 


36 


8 


-285 


-35 


160 


-146 


11 


85 


-368 


-0.1 


-16 


1392 


81 


36 












multi-coring 


urational, multi-reference 










3 


-221 


-52 


-52 


38 


-168 


-132 


-414 


-0.4 


-66 


1267 


-117 


25 


4 


-217 


-34 


166 


-106 


-24 


52 


-289 


-0.1 


-26 


1210 


52 


25 


5 


-246 


-37 


162 


-115 


-14 


59 


-306 


0.0 


-25 


1218 


71 


24 


6 


-242 


-38 


174 


-116 


-4 


52 


-298 


0.1 


-30 


1291 


78 


28 


7 


-266 


-42 


170 


-128 


6 


57 














final 


-272 


-38 


172 


-138 


10 


71 


-347 


-0.1 


-27 


1328 


86 


38 



The relativistic effects arc estimated from the differences 
between the non-relativistic CI and corresponding RCI 
results. 



In Table IV we compare the relativistic corrections on 



obtained with our calculations to the ones deduced 



from single configuration calculations (RCI-HF and DF- 
HF). It seems that correlation plays an important role 
in the estimation of these corrections, as could be ex- 
pected from an operator that measures the correlation 
between the momenta of the electrons. The convergence 
of the mono-reference approach with n is sufficient. How- 
ever, we see a large change between the mono- and multi- 
reference approaches. In Table [Vj we present the correc- 
tions for A-^ and B/Q that are both independent of the 
nuclear spin / and multipole moments 

Similarly to the non-relativistic calculations, we note 
that for neutral carbon, the impact of the 7 th and 8 th 
shells is not much affected by the choice of reference. We 
then estimate the final value as in the non-relativistic 
case. 



III. RESULTS AND COMPARISON TO 
EXPERIMENT 

A. Hyperflne Structures 

In this work, we focus on the isotopes 13 and 11 of 
carbon, respectively of nucleus spin 1/2 and 3/2. The 
n C nucleus decay into B by e + -emission with a half- 
lifetime of 20.4 minutes. Haberstroh et al. [JS] and Wol- 
ber et al. [44] performed experimental studies of the hy- 



TABLE VI. Aj- [kHz/^jv] and B/Q [kHz/barn] theoretical 
values for carbon 3 P, 1 D and C" 4 S° 2 D°. 





C 




cr 




state 


A— 

Mr 


B/Q state 


A— 


B/Q 






Non-relativistic 






2p 2 3 Pi 
2p 2 3 P 2 
2p 2 1 D 2 


2 296 
105 883 
156 317 


74 184 2p 3 4 Ss /2 
-148 367 2p 3 2 Dl /2 
283 206 2p 3 2 Dl /2 


9 394 
53 836 
107 317 




-35 455 
-50 650 


+ relativistic corrections 



2p 2 3 Pi 



2 024 



2p 2 3 P 2 106 055 
2p 2 X D% 156 327 



74 145 
-148 505 
283 276 



2p 3 4 S° 3/ 2 
2p 3 2 D° /2 
2/ 2 D° /2 



9 048 
53 809 -34 128 
107 403 -50 613 



pcrfine structures of the carbon ground state of 11 C and 
13 C, respectively. In the latter article, a magnetic dipole- 
moment of 11 C of —0.964(1) was deduced from the 
then available /x( 13 C) value. We update this estimation 
by using the modern /i( 13 C) value |47j combined with the 
two measured A( 3 P2) constants: 



M ( n C) 



A/( n C) In M 13 C) 



= -0.9642(2) fi N . 

(28) 

The error on this value is now dominated by the accuracy 
of the A( 3 P2) hyperfine constants measurements. 

As mentioned in Section |IIE| it is difficult to have a 
rigorous estimation of the uncertainty on the hyperfine 
parameters. We however advance a learned guess of their 
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TABLE VII. Comparison of our calculated hyperfine constants of C to other works. The experimental values are adjusted 
according to our analysis of the off-diagonal J J' interaction. All values are in MHz. 







13q 






i3 C - 






A^P) 


A^P) 


A 2 ( L D) 


A S/ 2CS°) 




A 6/2 r-D°) 


Original expQ 
This work 


2.838(17) 
2.84 


149.055(10) 
148.99 


219.61 


12.71 


75.59 


150.88 


Prev. worlj^] 


2.28 


148.1 











reliability. First, we see in Table |T] that the integrals ai, 
a s d and b q of C change less than 0.05% after the addition 
of the [~13f| correction. The a c parameter of the 3 P is 
only slightly more affected (~ 0.1%). These effects are 
representative of the accuracy of our results for neutral 
carbon. In the case of C - we face two additional limita- 
tions: the structure of the negative ion converges more 
slowly and we are limited in our expansions. Moreover, 
only the most troublesome contact term is responsible 
for the non-relativistic HFS of 2p 3 4 5°. For these rea- 
sons, and comparing the values of Table [TT] with results 
obtained with the active set |"10fc] , we must allow for rel- 
ative uncertainties on the HFS parameters roughly ten 
times larger for C~ than for C. 



are involved, one must diagonalize the matrix 



W(JJ';F) 
W(JJ';F) Ajj,E(LS F) 



(29) 



where Ajj,E(LS F) = E(LSJ'F) - E(LSJF) is dom- 
inated by the fine structure splitting and W{JJ'\ F) is 
governed by the off-diagonal hyperfine constants - here 
Ajj. 



IIBl The off- 



Bj,j-i and Bj t j_ 2 ( see Section 
diagonal electric quadrupole interaction is negligible and, 
at the non-relativistic level, we obtain for C( 3 P) 



IA^/hj = 50.47 MHz/^v , 
IA 2 ,i/hj = 62.71 MHz/ hn , 



In Table 



VI 



we present the non-relativistic A-^ and 



B/Q results calculated using the final values of a;, a se j, a c while for C ( 2 D°) we have 



and b q of Tables [T] and [II] In the same table, we add 
the relativistic corrections of Table |V] to those values. 
The A( 3 Pi) constant is the place of severe compensations 
between the orbital (a;) and spin-dipole (a s d) contribu- 
tions so that the uncertainties on those sum up to an 
error of the order of 10 2 kHz/fiN- The other nuclear- 
parameters-independent hyperfine constants of neutral 
carbon suffer of a non-relativistic uncertainty of about 
10 - 10 2 kHz//ijv These are larger than the fluctuations 
observed in Table [v] As far as C~ is concerned, on the 
one hand the A S° hyperfine structure is essentially due to 
the contact term, itself arising only from correlation ef- 
fects, and on the other hand, the 2 D° hyperfine constants 
are small but the achieved convergence of the calculations 
is less good. Therefore the relative non-relativistic uncer- 
tainties on C~ hyperfine structures are larger as they sum 
up to about 50— 100 kHz/^N (kHz/barn). We conclude 
that the reliability of all normalized hyperfine constants 
is of the order of 10 2 kHz/ hn (kHz/barn) with the ex- 
ception of the B( 4 S°) that is certainly negligible. 

O ur re sult s are compared with observations in Ta- 
bles VII and VIII for C and C~ respectively. We ob- 



IA 



5/2,3/2 



/Hi = 34.20 MHz//i 



N- 



(30) 
(31) 
(32) 



(33) 



The hyperfine interaction between states belonging to 
different terms is negligible. 

Wolber et al. [441 measured two hyperfine splittings in 
the 13 C 3 Pj multiplet, allowing the determination of the 
Ai and A 2 diagonal constants but not of the off-diagonal 
constants so that they had to deduce the contribution of 
the JJ'-interaction theoretically. The level shifts that 
they obtained from their computations are significantly 
higher than ours. However, the Aj constants that repro- 
duce the experimental hyperfine splittings when using 
our results for the JJ'-interaction, A 1 = 2.829(17) MHz 
and A 2 = 149.052(10) MHz, do not differ largely from 
the experimental constants presented in Table [VlTJ 

Haberstroh et al. [46 measured three hyperfine split- 
tings for n C 3 Pj, which is insufficient for determining 
all four Aj and Bj of this term. Hence they deduced the 
value of B\ from the relation B 2 /B\ = —2 which is only 
valid in the Hartree-Fock model. From Table |VT} we see 
that this formula holds very well at the non-relativistic 
level but that, including relativistic corrections, we have 



serve a good agreement with experiment, better than ex- 
pected from the above discussion. This represents a sig- 
nificant improvement compared to the theoretical study 
of Jonsson et al. [45j . 

The observed hyperfine splittings arise from the diago- 
nal hyperfine interaction, parametrized by the Aj and Bj 
constants, and, to higher order, from the non-diagonal 
(J J') interaction of states of same F. If only two levels 



B 2 /Bi 



-2.0029. 



(34) 



The effect of the refined B 2 /Bi ratio and JJ'-interactions 
cancel each other in the estimation of the diagonal hyper- 
fine constants so that the resulting Aj constants do no 
differ significantly from the experimental ones quoted in 
Table |VIII| For the electric quadrupole interaction, the 
accuracy of our results is such that we can safely update 
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TABLE VIII. Comparison of our calculated hyperfine constants of 11 C to other works. The experimental values are adjusted 
according to our analysis of the off-diagonal J J' interaction and of B ( 3 Pi) / B ( 3 P 2 ) . All values are in MHz. 









U C 










A^P) 


Pi( J P) 


A 2 ( J P) 


P 2 ( J P) 


A 2 ( L D) 


B 2 CD) 


Original expJ^J 


-1.308(24) 


2.475(14) 


-68.203(7) 


-4.949(28) 






This work 


-1.30 


2.474 


-68.17 


-4.955 


-100.49 


9.450 








ir c - 












B3/2CS ) 




B 3/2 ('D°) 


A 5/2 ( Z D°) 


B b/2 {'D°) 


This work 


5.82 


w 


34.59 


-1.139 


69.04 


-1.688 



the electric quadrupole moment of the 11 C nucleus with 
the formula 



Q("C) 



(5 2 ( U C)) 



( xp 



(B 2 /Q) 



(35) 



th 



Using the B2 constant of Haberstroh et al. [HO , we obtain 
a value of +0.03333(19)^(2)^ barns but if we use our 
theoretical parameters in the analysis of the observations, 
we obtain 



Q( n C) = +0.03336(19) ea: p(2) t/l barns 



(36) 



This value is used for estimating the theoretical Bj con- 
stants of this work presented in Table |VIII| The differ- 
ence between theory and experiment for the B2 constant 
follows directly from the fact that ( 36 1 includes the refine- 



ments of the theoretical parameters needed in the anal- 
ysis of the observed hyperfine splittings. 

Let us mention the previous calculations of the b q pa- 
rameter (we get b q = 0.6314 a^ 3 ): b q = 0.6325 a^ 3 gS], 
b q = 0.6319 cio 3 EH 



Our systematic procedure is not particularly efficient 
for predicting the negative ion binding energy. In par- 
ticular, for the 4 S° detachment threshold, the coupled- 
cluster approaches are much more impressive [151 116) . 
The recent value of Klopper et al. [16] indeed achieves a 
sub-meV (< 8 cm' 1 ) agreement with the experimental 
electron affinities for all first- and second-period atoms 
(H-Ne). A similar accuracy had already been achieved 
more than 10 years before by de Oliveira et al. [T5] for 
the second and third period p— block atoms. 

By trying various extrapolation schemes on our 
calculations, we explain up to ~ 20 cm -1 of the differ- 
ence between our calculation of the 4 S° binding energy 
and the experimental value (about 5 cm -1 for each n 
and I extrapolations, and about another 10 cm -1 for 
the extrapolation to a complete active set). Turning 
to the relativistic effects calculations, we see that the 
scalar contributions calculated with the CC methods give 
-21.54 cm" 1 [IS] and -22.83 cm" 1 Q2] while we obtain 
—37.95 cm -1 . The extrapolation being reliable to about 
a couple of tenth of percents and since the additional 



JJsing the experimental constant expected contributions are of the order of the cm -1 , we 

conclude that our BPCI relativistic corrections are still 
unbalanced. The problem of our relativistic corrections 
on the detachment thresholds is confirmed by the fact 
that, looking to Table [Till they are not well converged. 

Aside a possible unbalance in the relativistic effects 
estimation, our error is roughly proportional to the cor- 
relation contribution. We see that the differences be- 
tween the HF and experimental energy separations (see 
Table IX) are reproduced to ~ 0.1 — 0.7%, which is about 
the percentage of the C( 3 P) correlation energy we get. It 
means that p and n are good indicators of the percentage 
of the correlation effects included in a model. However, 
the uncertainty on our relativistic corrections and on the 
2 D° missing correlation is too large for an extrapolation 
based on this observation to be useful, e.g. for improv- 
ing the experimental determination of the position of the 



B 2 quoted in Table VIII Sundholm and Olsen gS] pro- 
posed Q( n C) = +0.03327(24) barns which would only 
tenuously agree with our estimation if the (B2{ 11 C)) exp 
value was to be improved. 

In the case of C~, the small 2 D° fine structure 
(1.75 cm -1 , see below), leads to JJ'-interaction shifts on 
the energy levels that are roughly 10 times larger than 
in the neutral atom ground term, i.e. of the order of 
0.1 MHz. 



B. Energy differences 



Table IX presents several calculated energy separations 
and compares them to other works. 

Our C _ term splitting is in very good agreement with 
experiment but, as will be seen below, this is partially 
accidental. 

Our results on the neutral atom ground configuration 
level spacings are systematically better than the ones of 
Froese-Fischer and Tachiev [49]. It indicates that, in this 
context, our relativistic corrections are reliable. For the 
3 P fine structure, we obtain as accurate results as recent 
fully relativistic calculations [50] . 



2 Dj levels. Indeed, a 0.5% uncertainty on our calculated 
correlation energies, which is no overestimation, reflects 
in corrections ranging from ~ 60 cm -1 in the case of the 
largest HF-experiment discrepancy, to about 7.5 cm -1 
for the 2 D° — 1 S threshold, i.e. of the same order of mag- 
nitude than the experimental uncertainty. 

We would like to stress another advantage of using the 
number of correlation layers n and p as parameters for 
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TABLE IX. Comparison of the theoretical and experimental energy level separations. The MCHF - CI calculations of the 
energy differences involving C~ 4 S° and 2 D° are obtained with p — 99.8% and 99.3%, respectivel y, w hile we take p = 99.95% 
for neutral carbon transitions energies. Relativistic corrections (+rel) are estimated from Table III and the +[13f| column 
corresponds to the final results of Tables |l| and |ll| All values are given in cm -1 . 



State 






This work 




Exp. 


Prev. 


Th. 






HF 


MCHF - CI 


+rel 


+ 1 13ZI 








C-( 4 5 3 ° /2 ) 
C-( 2 D° /2 ] 
C-( 2 D° /2 ] 




0.0 
14 541.11 


0.0 
9 936.73 


0.0 
9 921.01 
9 922.76 


0.0 
9 916.63 
9 918.39 


0.0 

9 913.5(82Q 






C( 3 P ) 
C( 3 Pi) 
C( 3 P 2 ) 

cCd 2 ) 
cCs ) 




0.0 
0.0 
0.0 
12 573.20 
30 508.75 


0.0 
0.0 
0.0 
10 193.42 
21 657.45 


0.0 
16.39 
43.31 
10 220.55 
21 691.79 


0.0 
16.39 
43.31 
10 213.76 
21 682.40 


0.0 

16.423 
43.41 b 
10 192.633 
21 648. 01 c 


0.0 

16.33 c 
43.033 
10 268. 23 c 
21 818. 60 c 


16.4 d 
43.33 


c~ 


C 






Detachment thresholds 






4 QO 

D 3/2 

4 QO 

°3/2 

4 QO 

D 3/2 
2 n o 
^3/2 

1J 3/2 ~ 

2 r>° 

U 3 2 


3 Po 
l D 
'So 
3 Po 
'D 2 
'So 


4 438.80 
17 012.00 
34 947.55 
-10 102.31 

2 470.89 
20 406.44 


10 200.50 

20 391.04 
31 853.18 

240.39 
10 440.27 

21 898.72 


10 132.96 

20 356.04 
31 825.28 

194.77 
10 420.99 

21 886.54 


10 141.49 

20 357.77 
31 824.42 

207.67 
10 427.10 

21 890.05 


10 179.68(16f] 

20 372.31(16/ 
31 827.69(16) f 

266.2(81) 
10 458.8(81) a 

21 914.2(81) a 


10 184.610 
43fl 


10 185.£D 



preserving the balance between systems having different 
numbers of electrons. Observing that our models con- 
verge toward the exact solution of the Schrodinger equa- 
tion, a larger number of electrons demands larger active 
sets, and that for a given orbital set, p is roughly pro- 
portional to the amount of correlation in the model, we 
affirm that the results obtained with increasing n and 
p will most often underestimate the photodetachmcnt 
thresholds. In other words, the detachment thresholds 
are valuable references for estimating the accuracy of the 
calculations since their behavior is monotone (as the level 
energies themselves). 



C. Mass isotope shifts 

Table [x] reports the results for the (A' = 13, A = 12) 
IS on various energy differences and compares them to 
previous works. 

We have seen that the S sms parameter is strongly cor- 
related to the energy, with negative and positive angular 
coefficients with respect to p and n, respectively. The 
non-monotonous behavior of the SMS forbids us to gen- 
erally conclude that any calculation similar to ours will 
result in upper or lower bounds to the differences in mass 
polarization expectation values. In our particular case 
however, Tables [T| and \U\ and Figure [2] show that the con- 
vergence in n is better achieved than in p (truncated to 
99.8% or 99.3%). Therefore, we likely overestimate S sms . 
Since we have globally AS sms > AE, the estimations of 
the IS on the detachment thresholds presented below, in 
particular on the 54, are probably overly negative. Fur- 



thermore, we estimate that our non-relativistic values of 
IS are reliable to about 2 lCT 1 nT 1 if the CT( 2 £>°) is 
involved and of the order of 10~ 2 m -1 if not. Hence 
we cannot explain the disagreement between the isotope 
shift on the 54 of Klopper et al. [16] and ours. This is not 
an isolated discrepancy since we can extract from their 
results an (18 — 16) IS on the 54 of oxygen of —11.7 m _1 
which is in disagreement with the experimental value of 
-7.4(18) m" 1 [58]. Godefroid and Froese-Fischer [23] 
obtained an IS of -5.73 m" 1 with a MCHF model, 
which is inside the experimental error bars (see also ref- 
erence 38J). The MCHF approach has also proven its 
usefulness for the calculation of IS on the 54 of heavier 
systems [251 125] . 

With regards to the difficulty to calculate the mass 
shifts, the discrepancy between the different calculations 
of the IS on the neutral atom term separation, of the 
order of 0.1 m _1 , is understandable. From the compari- 
son of our results and the ones of Kozlov et al. [50] , it is 
difficult to estimate an order of magnitude for the con- 
tribution of the relativistic effects. 

The neutral 3 P fine structure has been much more 
studied. It is known that the relativistic corrections to 
the mass shift operator are crucial when studying the iso- 
tope shift on the fine structure [S"5H6"2"] . Veseth [S3] and 
more recently Kozlov et al. |50j performed calculations of 
the relativistic mass shifts in the 3 P multiplet of carbon, 
the first by treating perturbatively the fine-structure and 
nucleus-mass dependent Hamiltonians up to the third or- 
der, the second by using an all-electron CI method on the 
Dirac-Breit Hamiltonian and calculating the expectation 
value of the relativistic MS operator valid to the second 
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TABLE X. Comparison of our isotope shifts (A' = 13, A — 12), in m~ , on various positive energy separations. The MCHF - 
CI calculations of the energy differences involving C~ 4 S° and 2 D° are obtained with p = 99.8% and 99.3%, respectively, while 
p — 99.95% for neutral carbon transitions. Relativistic corrections (+rel) are estimated from Table IV and the + |T3f| column 
corresponds to the final results of Tables |T] and [TTJ 









SMS 


NMS 


IS 




trans 




MCHF - CI 


+rcl 


+ |13Z| Exp. This work 


Other th. 


Exp 










IS on the C~ terms separations 






4 S'3/2 — 


^3/2 


-4.965 


-4.960 


-4.953 3.500 -1.454 
IS on the C~ detachment thresholds 






4 Q° 

°3/2 — 
2 n o 
U 3/2 


3 Po 
l Jh 


-10.656 
-8.865 


-10.629 
-8.875 


-10.630 3.592 -7.038 
-8.878 3.690 -5.185 

IS on the C terms separations 


-8.70 




3 P>- 
3 Po- 


l D 2 
'So 


-3.021 
-3.507 


-3.052 
-3.540 


-3.047 3.597 +0.550 
-3.536 7.640 +4.103 


+0.5053 

+4.672 b 





3 P1 - 3 P 2 



'Pi 



IS on the C( 3 P) fine structure 
-0.033 -0.033 0.010 -0.023 

+o.oi4f] 

-0.019 -0.019 0.006 -0.013 



+4.3743 



+0.020 b 
+0.015^ 
+0.009 b 
+0.010 f 



+0.0137(10 
+0.0180(43,, 
+0.0077(7£f 
+0.0057(83) 



order in aZ. We easily estimate the relativistic correc- 
tions to the specific mass shift operator for the transition 
3 Pi — 3 P 2 by comparing the equation (8) of Veseth to the 
Breit-Pauli fine-structure operator [53]. We obtain a cor- 
rection of +0.0375 m- 1 for the 13 - 12 P5( 3 Pi- 3 P 2 ) which, 
combined with our result of Table [Xj gives a total shift 
of +0.014 m _1 . This value is in good agreement with the 
observation and with Veseth's results. However, neither 
this observation nor the stability of the corrections of Ta- 
ble |IV] demonstrates that the scalar relativistic effects on 
the IS are reliable. 



D. Transition probabilities 

We study the Ml and E2 transition probabilities be- 
tween the LSJ levels of the ground configurations of C 
and C~. With the exception of 1 D 2 — 1 Sq and the tran- 
sitions between states of a same multiplet, the calculated 
Einstein coefficients are only non-zero thanks to the LS 

U b/2 



relativistic mixing. The Ml channel of the A S^j 2 
transition is itself only opened by LS mixing of corre- 
lation CSFs. The non-relativistic magnetic dipole and 
electric quadrupole interactions are computed using the 
BPCI wave functions based on the MRg 8 -I model de- 
scribed in the end of Section III Fl For the calculation 
of the Einstein A transition rates between states that 
are developed in non-orthogonal orbital sets, we use the 
biotr program that is part of the ATSP2K package [25] . 
The results of our calculations on the neutral atom 



are summarized in Table |XI| The convergence of the A 
coefficients with the active set is well achieved and the 
comparison with the results of Froese-Fischer [64 is fa- 
vorable. Still, we can point out that, for the weak E2 
transition rates of inter-term transitions ( 3 Po — 1 D 2 
and 3 Pi — 1 D 2 ), the relative change between the results 
of Froese-Fischer and ours is quite large. 

The results for the C~ are displayed in Table XII 
There is other value available in the literature. To fill 
this gap, we compare our transition probabilities with 
others for the first elements of its iso-electronic sequence, 
i.e. N I and O II. A priori, the omission of the relativis- 
tic corrections to the Ml transition operator JJ, 66 could 
be a serious limitation of our calculations. Indeed, Eiss- 
ner and Zeippen [67) showed that for transitions between 
terms of the 2p 3 configuration, in particular, the rela- 
tivistic corrections are of the same order of magnitude 
as the usual non r elativistic amplitude (see N I and O II 
data of Table pa!]). 

However, as can be seen from Table XII and in refer- 



ence [68], if the Ml channel becomes rapidly dominant 
with increasing Z along the nitrogen-like sequence, the 
E2 channel remains, for low Z, of the same order of mag- 
nitude as the Ml channel. Furthermore, in the case of the 
C~ system, the relativistic mixing of the 2 Dy 2 an( ^ 4 '-'3/2 
is even further suppressed by the diffuse nature of the 2 D° 
state. The Ml channel is then only a small correction to 
the total A in the CT ( 4 S° - 2 D°) transition, being of the 
same order of magnitude as the uncertainty with respect 
to the convergence with the active set (\8k~\ — \7i\). 
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TABLE XI. Einstein A coefficients, in s~\ for the 2p 2 intra-configuration Ml and E2 transitions of carbon calculated using 
the BPCI wave functions based on the MRgg-I models for the active space [lOfc] and \12k~] , compared with the calculations 
of [64]. Transitions vacuum wavelengths (A) are reported in Angstrom. 



This work 
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10 
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12 
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NIST 1651 
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A ki 
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A ki 
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3 Po 


- 3 Pi 
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6072(+3) 


8.033(- 


-8) 


6070(+3) 


8.041(- 


-8) 
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8.114(- 


-8) 
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3 Pi 


- 3 P 2 


Ml 


3703(+3) 


2.656(- 


-7) 


3701(43) 


2.660(- 


-7) 


3700(+3) 


2.662(- 


-7) 


3704(+3) 






E2 




3.529(- 


-15) 




3.536(- 


-15) 




3.633(- 


-15) 




3 Po 


- 3 P 2 
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2300(+3) 


1.696(- 


-14) 


2299(+3) 


1.700(- 


-14) 


2296(+3) 


1.754(- 


-14) 


2304(+3) 


3 P 2 


- 1 D 2 


Ml 


9805 


2.370(- 


-4) 


9814 


2.358(- 


-4) 


9735 


2.245(- 


-4) 


9853 


3 P1 




E2 




1.210(- 


-6) 




1.209(- 


-6) 




1.140(- 


-6) 






Ml 


9779 


7.513(- 


-5) 


9788 


7.504(- 


-5) 


9710 


7.544(- 


-5) 


9827 






E2 




1.050(- 


-7) 




1.037(- 


-7) 




1.576(- 


-7) 




3 Po 


- 'D 2 


E2 


9763 


7.789(- 


-8) 


9772 


7.836(- 


-8) 


9694 


6.242(- 


-8) 


9811 


3 P 2 


- 


E2 


4604 


2.138(- 


-5) 


4606 


2.130(- 


-5) 


4587 


2.250(- 


-5) 


4629 


3 P1 




Ml 


4598 


2.368(- 


-3) 


4600 


2.367(- 


-3) 


4581 


2.381(- 


-3) 


4623 


l D 2 


- 'So 


E2 


8679 


6.148(- 


-1) 


8681 


6.135(- 


-1) 








8730 



TABLE XII. Einstein A coefficients, in s _1 , for the intra-configuration Ml and E2 transitions of all C~ bound states, calculated 
using the BPCI wave functions based on the MRgg-I models for the active space [4/] to \8k] . We compare these results with the 
corresponding nitrogen and nitrogen-like oxygen A coefficients. Transitions vacuum wavelengths (A) are reported in Angstrom. 
The final set is obtained by re-normalizing the transition probabilities [~8fe] with these experimental ( 4 S° — 2 D°) and calculated 
( 2 2-5/2) energy differences. 
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n = 4 


8813 


1.177( 


-6) 


7.628(- 
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8815 


4.823(-8) 


1.118(- 


-6) 


-5716(+4) 


8.664(-ll) 


3.108(- 


-22) 


n = 5 


9356 


1.130( 


-6) 


8.915(- 


-7) 


9357 


4.197(-8) 


1.297(- 


-6) 


-1245(+5) 


8.415(-12) 


3.912(- 


-23) 


n = 6 


9590 


1.126( 


-6) 


1.019(- 


-6) 


9590 


3.999(-8) 


1.470(- 


-6) 


-6509(+5) 


5.893(-14) 


2.141(- 


-26) 


n = 7 


9801 


1.123( 


-6) 


1.125(- 
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9800 


3.802(-8) 


1.605(- 
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1519(+5) 


3.092(-12) 


4.209(- 


-23) 


n = 8 
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1.080( 
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1.688(- 
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-22) 
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1.046(-6) 
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1.414( 
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3.382(- 
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1.241(-7) 
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1.45(- 
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2.13(- 
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7.58(-6) 
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5) 
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1.30(-7) 






Z87 c 




1.58(- 


4) 








2.00(-6) 








1.30(-7) 







The inversion of the fine structure splitting of the 2 D° 
term in C~ with increasing active set reflects the diffi- 
culty to calculate this quantity for half filled shell sys- 
tems, particularly in highly correlated systems. From 
this regard, it is also interesting to note that this fine 
structure is "normal" , i. e. not inverted as in the heavier 
iso-electronic systems. In other words, it tends to behave 
like a less-than-half-fillcd-shell system. 

The experimental ( 4 Sy 2 — 2 D°j) transition is at 

10087(9) A, and the 2 D° calculated fine structure that 
we recommend (see Table IX I is 4-5703 10 4 A. Ta- 



ble |XII| presents the final A coefficients calculated with 
the [8ft] active set and renormalized by the experimental 
( 4 S° — 2 Dj) energy separation, and our t heoretical value 
for the 2 D° fine structure (see Table IX I. 



IV. CONCLUSION 

We performed large scale MCHF-CI calculations of the 
energy levels belonging to the lowest configuration of neu- 
tral carbon and all bound states of C~, including the fine 



14 



structures, hyperfine structures and isotope shifts. In 
addition, we calculated all Ml and El transition rates 
between the studied LSJ states. 

The overall precision of the non-relativistic expectation 
values is estimated to be about 0.3 — 0.8%. However, 
this imprecision on the total energy and S sms leads to 
a larger uncertainty on the differential effects. To our 
knowledge, the obtained non-relativistic energies are the 
most accurate ah initio values to date. We obtained an 
extrapolated energy for the C( 3 P) state of -37.84465 E h 
in good agreement with the semi-empirical "exact" value 
of -37.8450 E h 02]. 

We conducted a careful study of the relativistic cor- 
rections deduced from the comparison of relativistic CI 
calculations with the corresponding non-relativistic cal- 
culations. Even if the so-deduced corrections do permit 
a relevant theory versus experiment comparison and that 
the correlation effects still dominate our uncertainties in 
many cases, an estimation of the relativistic effects on a 
firmer basis should be performed in more accurate stud- 
ies. We also note that the C~ negative ion is very little 
affected by relativity. 

We check the experimental hyperfine structure of the 
neutral carbon by replacing the theoretical parameters 
used in the original papers [44j [46] by our values. The 
resulting hyperfine constants and 11 C nuclear magnetic 
moments are in good agreement with previous experi- 
mental and theoretical studies. 

As far as the transitions probabilities calculations are 



concerned, we find a good agreement of our neutral car- 
bon A coefficients with the ones of the literature. For 
the C - intra-configuration transitions, we expect the rel- 
ativistic corrections to the Ml operator to be less impor- 
tant than in higher Z iso-electronic systems. Once more, 
the missing correlation effects are equally limiting. 

We find that the parametrization of the model in terms 
of the number of correlation layers (~ n) and percent- 
age of the wave function accounted by the MR in sub- 
sequent CI calculations (= p), provides useful tools for 
including a fixed percentage of the total correlation ef- 
fects. It allows to establish lower bounds on detach- 
ment thresholds and, for calculations that are sufficiently 
converged with respect to n, upper bounds on AS sms 
(lower bound on S sms ). 
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